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Abstract 

A relationship between the time scales of quantum coherence loss and short-time 
solvent response for a solute/bath system is derived for a Gaussian wave packet approx- 
imation for the bath. Decoherence and solvent response times are shown to be directly 
proportional to each other, with the proportionality coefficient given by the ratio of 
the thermal energy fluctuations to the fluctuations in the system-bath coupling. The 
relationship allows the prediction of decoherence times for condensed phase chemical 
systems from well developed experimental methods. 

PACS numbers: 31.70.Dk, 42.50.Lc, 78.47.+p, 82.20.Wt 

Quantum processes in condensed phases are often studied by focusing on a small subset 
of degrees of freedom and treating the rest as a bath. The subsystem of interest may com- 
prise a single molecule, a molecule plus its nearest surroundings, or even a single vibrational 
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mode within a molecule. The remaining degrees of freedom form the bath. In the presence 
of system-bath interactions, the subsystem's wave function evolves into a superposition of 
quantum states. Due to an enormous density of states in a macroscopic environment, small 
differences in the system-bath coupling lead to rapid divergence between bath evolutions 
corresponding to different states of the subsystem. The reduced density matrix of the micro- 
scopic subsystem, obtained from the total density matrix by integrating over bath degrees of 
freedom soon b ec 0m es diagonal Quanta states decohere. The decohere rate is deter- 
mined by the sensitivity of bath evolutions to the quantum state of the subsystem. Notably, 
this sensitivity also determines bath response to a perturbation within the subsystem. In the 
context of the condensed phase chemical physics, the rate of the solvent bath rearrangements 
following a perturbation of the solute subsystem is described by a well developed solvent re- 
sponse theory.! - ! In this Letter, we establish for the first time a quantitative relationship 
between quantum coherence loss and the short time solvent response. 

The new relationship is important because a link between solvation dynamics and deco- 
herence is capable of providing valuable insights into both phenomena. Modern techniques 
of the solvent response theory employing the concepts of dielectric and mechanical relax- 
ation,!'! and instantaneous normal mo des9 then become transferable to the description of 
decoherence. The theory of short-time solvation can benefit from the recent theoretical ideas 
on quantum Brownian motionl'l and quantum measurement ,Bef0 where the notion of deco- 
herence appeared first. Most importantly, currently available experimental means to measure 
quantum coherence lossEMl and short time solvent response!@0 in a non-equilibrium sys- 
tem can be combined. A relationship between decoherence and solvation time scales further 
provides a tool to deconvolute the contribution of each effect on observed dynamics. The 
results presented below form a basis for evaluation of quantum decoherence times in various 
solute-solvent systems based on the extensive solvent response data accessible from both 
experimental measurements and adiabatic molecular dynamics simulations. 

Following Ref. ITS], we consider the direct product of the system S and bath B Hilbert 
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spaces. For simplicity, the system space is assumed to be two dimensional. We consider 
two orthogonal states <j)f and </>f of the system and a system-bath interaction that induces 
quantum transitions in the combined system, with the bath state responding to that of the 
system: 

<t>l®<t>* ^<t>l®<t>l, « = 1,2, (1) 

where <pQ is the initial state of the bath. Transitions from an arbitrary initial system state 
Ci4>f + c 20f are then described in terms of the reduced density matrix 
( Icil 2 c lC * 2 \ ( | Cl | 2 c lC ^f|0f) \ 

\ c \c 2 ) ^ \ c \c 2 {^m | C2 | 2 ^ 

Decoherence is defined^ as decay of the non-diagonal matrix elements, which, for the reduced 
density matrix of Eq. (||]), is clearly determined by the decay of the inner product of the bath 
states (0f Initially, the bath wave functions coincide: <pf — (j) 2 = 4>q. Later on, 

the bath wave functions correlated with the different states of the system diverge, and the 
overlap integral decreases. It is not the decay of the non-diagonal matrix elements per se 
that is most important from the practical point of view, but rather it is the associated 
slowing down of quantum transitions,^^ known as the quantum Zeno effect in the limit of 
infinitely fast decoherence.0 The life time of the quantum state in the presence of a bath 
varies inversely with the decoherence time. This result follows from the Fermi golden rule 
in the spin-boson model, where the bath is treated as a set of harmonic oscillators,0'El or in 
the frozen Gaussian formulation,&H where the bath wave function is approximated by a set 
of Gaussian wave packets.il Within this Gaussian wavepacket approximation, the average 
decay of the overlap integral and corresponding non-diagonal matrix elements in Eq. (0) is 



described by the decoherence function given by Eq. (39) of Ref. |2 



D(t) = (0fl0f) = ex P [-($: i ^)t 2 ], (3) 

n n 

where AF n = F\ n — F 2n is the expectation value of the difference in the quantum forces 
experienced by the nth bath degree of freedom, and the angular brackets indicate thermal 
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averaging. We note that here we do not consider the dynamical effect of the bath on the 
subsystem energy eigenvalues. In the presence of bath induced fluctuations in these values, 
a standard pure dephasing contribution (see Refs. |22l, |24]) to the decay of the off-diagonal 
elements of the reduced density matrix can potentially contribute as well. In the adiabatic 
representation the forces are given by the Hellmann-Feynman theorem 



V n #|0f), a = 1,2. (4) 



— 1/2 

In the low temperature regime the width a n of the Gaussian wave packet equals the 
width of the coherent state of a corresponding harmonic oscillator, i.e., a n = m n uj n /h. For 
higher temperatures, the width of the wave packet incorporates quantum thermal ensemble 
averaging. The thermal width is analytic for harmonic baths 

m ^tanhfeV (5) 



h \ 2k b T 



For arbitrary baths, the width can be defined via the thermal de Broglie wave length Xb 
(2iih 2 / inksT) 1 / 2 . An alternative expression for the thermal width is derived in Ref. ^0 



Sec. IV by comparing the exact and Gaussian wave packet results for the transition rate in 
a double well system 

Tfl n U/ n . , . 

■An, (6) 



A 



h 



coth 



2k B T ) hjj n 

This expression reduces to the coherent state width in the low temperature case, and gives 

in the high temperature limit. The last formula is particularly useful, since it yields a width 
which is independent of the frequency; Eq. (|7|) is designed^ for use in molecular dynamics 
simulations, where thermal averaging over bath states is performed classically. 
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Turning to solvation dynamics, the response of the solvent bath to a quantum transition 
within the solute subsystem is quantified by the normalized correlation function C of the 
energy gap U between the two quantum states.BB The fluctuation-dissipation theorem 
relates the non-equilibrium solvent response to the regression of fluctuations 5U of the gap 
U in equilibrium 

(5U(t) ■ sum 

C ® - m?) ■ ( 8 ) 

The short time solvation dynamics that is of relevance in the present discussion depends 
solely on the change in the solute-solvent coupling due to the quantum transition. The 
microscopic short time expression for C(t) has been obtained in Ref. [25| [Eqs. (2.18), (2.19)] 
by expanding C(t) in a set of independent modes and in time, yielding 



C(t) 



exp 



knT 




2{{5Uf) 

where U' n is the derivative of U with respect to the nth mass-weighted solvent coordinate 



(9) 



UL 



m. 



-1/2 



dU 
dx 



m. 



- 1/2 AF n 



(10) 



The decoherence t d and Gaussian solvation r 9 time scales are given by the variances of 
the decoherence D(t) and solvent response C(t) functions of Eqs. (|]) and (0), respectively. 
The structure of the equations is clearly similar, and by comparison we obtain 

\ \ f ■ ~ i ill \ 



T, 



9 J 



2k B T En Ai? nK 



((^) 2 )E„A^(a n ^)-i 
With the high temperature limit expression for the width [Eq. ([j])] the formula simplifies to 

12 1fl /7. rn\2 



Tn 



12(k B T) 2 _ Qk B T 
((SU) 2 ) ~~' 



(12) 



Here, 2A is the Stokes' shift, the difference between the equilibrium absorption and emis- 
sion maxima. Within a linear response regime, A is relatedll! to the fluctuations in the 
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quantum energy gap by the fluctuation-dissipation theorem: ((5U) 2 } = 2Xk B T. Eq. 
establishes direct proportionality between the short time evolution of decoherence and of 
solvent response in the high temperature limit for the bath. The proportionality coefficient 
is determined by the ratio of the thermal energy fluctuations (ksT) 2 to the fluctuations in 
the system-bath coupling ((5U) 2 ). The letter can be expressed in terms of the Stokes' shift 
in the linear response regime. We note that the dependences are sensible: High temperatures 
accelerate solvation dynamics, large equilibrium fluctuations in the coupling, as well as large 
Stokes' shifts, are indicative of fast decoherence. 

We note that the relationship between the decoherence and solvation time scales pre- 
sented above using the Gaussian wave packet approximation for the bath wave function 
also necessarily pertains to the spin-boson model. The spin-boson Hamiltonian describes a 
two-level system linearly coupled to a harmonic bathEi 



H S B = ~-hAa x + -Uq<T z + ^ \\ m nU n X 2 n + ^T) + \ < & (T * 



(13) 



where, Uq and q$ are the energy and coordinate displacements between the pair of potential 
minima, A is the intrinsic coupling between the two quantum states, c n is the system-bath 
coupling constant, and a z , a x are the Pauli matrices. The terms containing a z describe the 
energy gap 

U SB = U + q ^2 

(14) 

n 

According to Eqs. (||), (10), the short-time solvent response function of the spin-boson model 



is 



n 



The decoherence function can be extracted from the Fermi golden rule result for the spin- 
boson problem [Eqs. (3.35), (3.36) and (3.2) of Ref. 



D SB = exp 



(q 2 /Tih)Q 2 (t) , (16) 
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with 

Q 2 (t) = [°° J ^ )il ~ ' C ° SUjt) aoih(hu/2k B T)du, (17) 
Jo 

and 

The short-time expansion of Q 2 {t) gives 
D SB = exp 

l 4 * — * r 

n 

1 ^0 C n + 2 



- Y -^5- coth(^ n /2A; B T)t 2 l (19) 
4 ^— ' m n nw n -I 



exp 



4 ^ a„fr 2 



(20) 



with a n as in Eq. ([]). The expression in Eq. Q20| ) is a specific case of Eq. (g) with AF n = 
—qoC n . Comparison of Eqs. flT5| ) and (EO) leads to the relationship given in Eq. ( |TTj ) between 
the decoherence and Gaussian solvation times. 



We can directly test Eq. (12) for the case of relaxation following transition from the first 
excited to the ground state of the hydrated electron, whose time dependent properties are 
well studied theoretically (see Refs. ^-|22],|26|,p7] and references therein). Based on a molecu- 
lar dynamics trajectory for the first excited state of the hydrated electron,^ the equilibrium 
energy gap U is 0.56 eV and the fluctuation in the energy gap ((SU) 2 ) 1/2 is 0.21 eV, which 
corresponds to a 1.7 eV Stokes' shift (2A) at room temperature. This Stokes' shift corre- 
sponds closely to that found in non-adiabatic simulations.il The short-time component of 
the solvent response function C(t) of Eq. (§) is found to be characterized by a 10.6 fs Gaus- 
sian time scale, r g . From these data, the decoherence time t d of the first excited state of the 



hydrated electron calculated via Eq. fll2|) is 4.5 fs. This estimate falls within the previously 
reported range of 2.7-5.1 fs,&il with the value of 5.1 fs obtained in the high temperature 
approximation for Eq. (^). 

The properties of the solute and the nature of the quantum transition define the difference 
in the solute-solvent potential for the initial and final states and, therefore, determine the 
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magnitude of the Stokes' shift in the relationship flT2|) between the decoherence and solvation 
times. Intrinsic solvent properties also, to some extent, affect the magnitude of the Stokes' 
shift. However, the major solvent influence on the duration of quantum coherence is due 
to solvent's ability to respond to a perturbation in the solute, i.e. as reflected in the rate 
of solvent response, r g . Based on the success above for the hydrated electron, one can 
address other systems, which will demonstrate the variability of decoherence times. For 
a styryl dye in methanol,il experiment yields an estimated r g of 40 fs and a Stokes' shift 
of 115 nm yielding a decoherence time of 6.8 fs. Acetonitrile [CH3-CN], the next solvent 
in a logical series exhibits a 100 fs experimental short time solvent response.^ Simulation 
of the electronic transition of the betaine-30 molecule in acetonitrile is characterized by a 
similar 91 fs short-time solvation. The value for r g , together with the 0.16 eV electronic 
energy gap fluctuation evaluated along the ground state trajectory ,0 leads via Eq. fli~2|) to a 
substantially longer 49 fs decoherence time. Compared to the protic solvents, acetonitrile is 
much less effective in destroying quantum coherence. 

In summary, we have presented an analytical relationship between the time scale for 
quantum decoherence and that governing the short time response of solvent to a pertur- 
bation in the solvent-solute coupling. The proportionality constant relating these requires 
knowledge of only the Stokes' shift associated with the change in solute state. The expres- 
sion successfully reproduces results obtained directly by other routes, allows the prediction 
of decoherence times for other solution systems. It is expected that since the required in- 
put data is becoming readily accessible experimentally for even the most rapidly responding 
condensed phase environments, the derived relation will be very valuable to advancing the 
study of both decoherence and condensed phase chemical dynamics more generally. 
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